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ABSTRACT 

We present the results of a multi-site photometric campaign carried out in 2004-2008 for the Algol-type eclipsing binary 
system CT Her, the primary component of which shows S Scuti-type oscillations. Our data consist of differential light 
curves collected in the filters B and V which have been analysed using the method of Wilson-Devinney (Phoebe). 
After identification of an adequate binary model and removal of the best-matching light curve solution, we performed a 
Fourier analysis of the residual B and V light curves to investigate the pulsational behaviour. We confirm the presence 
of rapid pulsations with a main period of 27.2 min. Up to eight significant frequencies with semi-amplitudes in the 
range 3 to 1 mmag were detected, all of which lie in the frequency range 43.5-53.5 . This result is independent 
from the choice of the primary's effective temperature (8200 or 8700 K) since the light curve models of the binary 
are very similar in both cases. This is yet another case of a complex frequency spectrum observed for an accreting S 
Scuti-type star (after Y Cam). In addition, we demonstrate that the amplitudes of several pulsation frequencies show 
evidence of variability on time scales as short as 1-2 years, perhaps even less. Moreover, our analysis takes into account 
some recently acquired spectra, from which we obtained the corresponding radial velocities for the years 2007-2009. 
Investigation of the 0-C diagram shows that further monitoring of the epochs of eclipse minima of CT Her will cast a 
new light on the evolution of its orbital period. 

Key words, stars: binaries: eclipsing - stars: oscillations - stars: fundamental parameters - stars: individual: CT Her 

1. Introduction ily in low overtone radial and non-radial acoustic modes 

, . . , 1 1 • 1 • with periods between 30 min and 6 hr. Some may also pul- 

H . Asteroseismology aims to understand the pulsation physics ^^^^ ^^^^-^^ ^^^^^ (|Kurtz 2000^ . Their excitation mech- 

C?; m order to probe the interior parts of all kinds of stars. ^^^^^^ ^ ^ «:-mechanism active in the partial ionisa- 

Among the classical pulsators, S Scuti stars are rather com- ^^^^ ^^^^^ jj^ j jj^ jj^ well-understood. However, 

mon and are located at the intersection of the Cepheid m- ^^-^^ ^^^^^ unknown amplitude limiting and mode selec- 

stabihty strip and the mam sequence. They pulsate primar- mechanisms are operating in these pulsators. Only a 

7^ ; ^ ~~ ~~ T fraction of the theoretically predicted modes are observed, 

bend offprint requests to: P. Lampens , . , i, • /, n r , • ,i i 

* „ , , , , . , , ,,,,,,, , , . which results m too many tree parameters m the pulsa- 

Based on photometric data collected at the observatories ,. i i t, • • i i • i- i , n , i i i 

listed in Table 1 and spectra acquired at the NAO, Bulgaria, and tion models. It is indeed primordial to collect knowledge on 

Calar Alto Observatory Spain. The Skinakas Observatory is a t^e fundamental stellar properties of the pulsator mdepen- 

collaborative project of the University of Crete, the Foundation dently. 
for Research and Technology - Hellas, and the Max-Planck- 

Institut fiir Extraterrestrische Physik. One way is to study the pulsating components of binary 

** Tables [3l H [9] and [12] are only available in elec- (multiple) systems in great detail, since binary (multiple) 

tronic form via anonymous ftp at cdsarc.u-strasbg.fr or at systems with well-characterized components supply addi- 

[http://csdweb.u-strasbg.fr/cgi-bin/qcat?/A-|-A/XXX/YYY( tional constraints for a more reliable modelling. The posi- 
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tions in the H-R diagram and therefore the components' 
evolutionary statuses can be more accurately determined 
than in the case of single stars, e.g. for those 6 Scuti stars 
that are at the (very) end of their H-core burning phase. 
The complexity of having to deal with additional compo- 
nents does generally not weigh against the scientific return, 
although observations spread over even longer time-scales 
may be necessary to disentangle both phenomena. 

About seventy percent of all stars of the Solar neigh- 
bourhood are members of binary or multiple systems 
(67% for G-M stars, [Mayor et al. 2001[ 75% for 0-B 
stars, IMason et al. 200ip . Yet, these facts are usually ig- 
nored in the study of stellar pulsation. There is, how- 
ever, strong evidence that duplicity affects the pulsa- 
tion properties in specific cases (e.g. the eccentric bina- 
ries HP 177863 (|De Cat fc Aerts 2002)) and HD 209295 
([Handler et al. 2002)) ). This is also predicted from a the- 
oretical point-of-view dWitte fc Savonije 1999[ ). It is thus 
essential to understand the possible link(s) between bina- 
rity and pulsation(s) and to observe whether or not - under 
given circumstances - the internal structure and the pulsat- 
ing properties of such stars might be different from those 
of the single pulsators ( [Lampens 2006[ ) . 

Detached and semi-detached eclipsing binaries (EEs) 
are particularly powerful tools in astrophysics: their accu- 
rate observation enables to derive the fundamental proper- 
ties (masses, radii, luminosities) of each component. EEs 
also provide the component's effective temperatures, as 
well as the distance if stellar atmosphere models are used 
([Maceroni 2006^ . 

We present a detailed photometric study of the Algol- 
type binary CT Her, an eclipsing binary of mag 11- 
12 and spectral type A3V-|-[G3IV] with an orbital pe- 
riod of 1.7863748 days ([Samus et al. 2009^ . Radial ve- 
locity measurements have also been collected, but due 
to its faintness, high-quality spectra are not easily ac- 
quired. CT Her belongs to the group of oscillating 
Algol-type (oEA) stars comprising w 35 known members. 
Up to now, a few members of this group only have 
been investigated thoroughly: RZ Gas has been observed 
photometrically as well as spectroscopically over almost 
a decade ([Lehmann fc Mkrtichian 20041 Rodriguez et al. 
IMHl [Soydugan erar2006l iLehmann fc Mkrtichian 2008[) 
while a very detailed photometric study of Y Cam 
was recently finished ( [Rodriguez et al. 2010[ ). TW Dra 
is another case for which a large spectroscopic ob- 
servational effort was carried out ([Lehmann et al. 20091 
ITkachenko et al. 2010[) . Its (0-C) diagram shows a slow 
orbital evolution ([Kreiner et al. 20"0T)l . The primary com- 
ponent displays oscillations of type S Scuti with a (main) 
pulsation period of 0.46 hr (« 28 min) and a total amplitude 
of about 0.02 mag (Kim et al. 2004a). Among the currently 
known oEA stars, it has one of the shortest orbital periods 
and the highest ratio Porb/-Ppuis (about 95). 

The oEA stars are former secondaries of evolved, semi- 
detached eclipsing binaries which are (still) undergoing 
mass transfer and form a recently detected subclass of 
pulsators close to the main sequence (Mkrtichian et al. 
120021 (2004)) . Searches for new oscillating Algols have been 
performed by Kim et al. ([2 004b, 20 01) and Mkrti chian 
et al. (p nrim [MlH f. [P igulski fc Michalska (2007) ' and 
Michalska fc Pigulski (2008 ) looked for them in the ASAS- 
3 and the UCLE-ll public databases. A search using the 
NSVS database has also been on-going at the Institute 



of Astronomy of the Bulgarian Academy of Sciences 
(Dimitrov et al. I2008al I2008bl I2009al I2009bp . Their 
general characteristics and pulsational properties have 
been summarized by Mkrtichian et al. (|2005p . Such oEA 
stars are, indeed, excellent laboratories for investigating 
the effects of mass accretion events as well as of tides 
unto the pulsation properties. Changes of these properties 
(amplitudes, modes and/or phases) due to mass accretion 
episodes have been observed (e.g. in RZ Cas where strong 
modal amplitude variations followed an abrupt change of 
its orbital period; Rodriguez et al. 120041 Mkrtichian et 
al. '2005). Their pulsational frequencies could be tidally 
split (as in KW Aur, IFitch 1976) and/or there may be 
some coupling between the pulsation and the orbital 
frequency due to some resonance mechanism (Mkrtichian 
et al. I2005p . Their very different evolutionary history is a 
challenge for stellar evolution modelling. For example, in 
order to reproduce the gainers of RZ Cas, KO Aql and 
S Equ through conservative binary evolution, large initial 
mass ratios (typically > 3) would be necessary. Eut the 
corresponding high mass loss rates in the beginning of the 
Roche Lobe overflow stage result in radii larger than the 
Roche radii. Hence, a non-conservative approach is needed, 
in which processes are considered that enhance the period 
without losing too much mass ([De Greve et al. 2009p . 
Such studies have not yet been attempted. Detailed and 
multi-year studies of more oEA stars are essential in order 
to provide solid grounds for their (future) asteroseismic 
modelling. 



2. Observations and data reduction 

2.1. The multi-site campaign 

High-precision light curves of CT Her (GSC 01509-1142, 
V = 11.347, B - V = 0.203) were collected in the frame- 
work of a 5-yr long multi-site campaign set up with the 
purpose to study its pulsational behaviour. A logbook of 
the observations is given in Table [TJ We observed in a dif- 
ferential mode from late spring till late summer of the years 
2004-2008. The original observations which led to its dis- 
covery as an oEA star (Kim et al. [2004a) were also included. 
The following comparison stars were used: C1=GSC 01509- 
1140 {V = 10.951,5 - V = 1.414); G2=GSC 01509-1052 
{V = 11.405,5 - V = 0.779); G2'=GSC 01509-1130 {V = 
12.02); G3=GSC 01509-0901 {V = 12.20,5- V = 0.78) 
and C4=GSC 01509-1090 {V = 10.75, B-V = 0.48). 

The characteristics of the various CCD cameras are the 
following ones: 

- at the Skinakas Observatory of the University of Crete, 
Greece. The camera is a Photometries 1024 x 1024 with a 
SITe SI003E chip of grade 1 and a pixel size of 24 fim cor- 
responding to a scale of 0.5" on the sky. The field-of-view 
is 8.5' X 8.5' ( [Papadakis et al. 2003p . 

- at the Sobaeksan Optical Astronomy Observatory of 
the Korea Astronomy and Space Science Institute (KASI), 
South-Korea. The camera has a 2048 x 2048 SITe chip with 
a pixel size of 24 /^m corresponding to a scale of 0.60" on 
the sky. The field-of-view is 20.5' x 20.5'. 

- at the Observatory of Mt. Lemmon, Arizona, operated by 
KASI. The camera has a 2084 x 2084 Kodak chip with a 
pixel size of 24 fim corresponding to a scale of 0.64" on the 
sky. The field-of-view is 22.2' x 22.2' large. 
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- at the Observatorio de MonegrillcQ, North of Spain. The 
camera is a SX Starlight CCD with a Sony ICX027BL 
chip (cooled to about -25°C) and a pixel size of 12.7 fi x 
16.6 n corresponding to 1.80" x 1.38". The field-of-view 
covered a sky region of 11.5' x 7.7'. The reduction was 
done using a software package called LAIA (Laboratory 
for Astronomical Image Analysis) developed by Joan A. 
Cancfl. 

- at the Beersel Hils Observatory, Belgium. The camera is 
a SBIG STlOXMe with a chip of grade 1 and a pixel size 
of 6.8 /xm corresponding to a scale of 1.43" on the sky (in 
2x2 binning). The field-of-view is 17.5' x 26' on the sky. 
The reduction was done using the Mira AP (v. 7) packago 

- at the Observatorio Sierra Nevada, South of Spain. The 
camera has a 2kx2k chip with a pixel size of 13.5 fim cor- 
responding to a scale of 0.23" on the sky (2x2 binning). 
The field-of-view is 7.92' x 7.92' large. 

- at the National Astronomical Observatory (NAO) 
Rozhen, Bulgaria. The CCD is a VersArray 1330B with 
a 1340 X 1300 E2V CCD36-40 chip of grade 2 and a pixel 
size of 20 /im corresponding to a scale of 0.258" on the sky. 
The field-of-view is 5.76' x 5.59' large. 

Photoelectric photometric data on CT Her were also 
acquired in the four Stromgren passbands uvby at the 
Observatorio Sierra Nevada between 2008 March, 31 and 
2008 April, 13. The comparison stars were Kl = HD 145122 
{V = 6.13), K2 = HD 146101 {V = 8.24) and K3 = HD 
145549 {V = 8.16)). The observations were carried out 
in the sequence Sky, Kl, K2, CT Her, AO Ser, Kl, K3, 
CT Her, AO Ser. . . because both oEA stars are close enough 
to one another on the sky. Unluckily, the orbital period was 
not fully covered during the additional observations. 

2.2. Observational technique 

CT Her was observed using the standard (Johnson) fil- 
ters B (mostly) and V (less frequently). At the Beersel 
Hills Observatory (BHO), we follow the specifications for 
the filters of Bessell (1995). CI is the principal comparison 
star which was commonly used at every observatory. We 
adopted the filter B as the main filter because of its higher 
signal-to-noise ratio and the larger amplitudes of pulsation 
expected. Typical exposure times for the B filter were set 
between 15 s (e.g. Sierra Nevada) and 60 s (e.g. Skinakas). 
Typical exposures for the V filter were 20-30 s (e.g. BHO 
and Skinakas) . All observers followed the standard calibra- 
tion procedure: a set of biases (resp. darks) was taken regu- 
larly during each night and a set of 5 to 6 fiat-fields per filter 
was obtained during evening and/or morning twilights. 

As an example, we describe the full reduction procedure 
used for images collected at the Skinakas Observatory. All 
the primary reduction steps were performed using stan- 
dard ESO-MIDAS routines. The frames were processed 
as follows: subtraction of the residual bias pattern using 
a median master bias frame, flat-fielding using a median 
master fiat-field frame, and median cosmic ray cleaning. 
Since the field is not crowded, the technique of aperture 
photometry was applied to extract the differential magni- 
tudes. The fixed aperture photometry was performed using 

^ http:/ /www. astrogea.org/jvidal/index. html 
^ http:/ /www. astrogea.org/soft/laia/laia.htm 
^ The Mira software is a registered trademark of Mirametrics, 
Inc., (http: / /www. mirametrics.com/index.htmi 



DAOPHOT (|Stetson 1987|) . CT Her and the comparison 
stars CI, C2' and C3 were measured using an aperture size 
as close as possible to the value providing the highest signal- 
to- noise ratio dStrigachev 2009[ ). The data consist of differ- 
ential photometry of the variable star in the sense (CT Her 
- CI) and of the check stars in the sense (C2'- CI) and (C3 
-CI). 

CT Her was also observed in the BVR filters together 
with standard fields (|Landolt 1992p . We performed all- 
sky photometry up to an airmass of 2 at the Skinakas 
Observatory on Aug. 2005, 1. The coordinates, calibrated 
magnitudes and colours of the target and comparison stars 
are listed in Table [2] 

3. Photometric data - time series 

The influence of interstellar and atmospheric extinction 
may be crucial when constructing light curves of eclips- 
ing binaries (jPrsa fc Zwitter 2005a ). This is generally the 
case when matching data from different sites and cam- 
paigns but it may also be valid when using data from a 
single site. Observations from different nights and/or sites 
do usually not match - their mean level may be shifted, 
and there may also a correlation with airmass during some 
nights. In this case, differential corrections associated with 
the second-order coefficient of the atmospheric extinction 
(k" ) were needed since the colours of the comparison star 
(CI) and the variable star are not the same (see Table 2). 
These corrections were mainly needed for the data acquired 
in the filter B because of the dependence in wavelength and 
the greater inhomogeneity of the i?-data sets collected over 
a longer period and at a larger number of sites equipped 
differently. 

To correct for this infiuence, we first computed the mag- 
nitude of CT Her based on the magnitude of CI in Tabled 
and the differential values (CT Her - CI). Next, we adopted 
a preliminary model light curve based on known and fit- 
ted parameters for the V light curve obtained during the 
2005 & 2008 runs at the Observatory of Monegrillo, Spain 
(cf. Table [1]). This was possible because the V light curve is 
less dependent on these corrections. In general, a clear de- 
pendence of the residuals in magnitude difference with the 
airmass per night was found, which we modelled using lin- 
ear regression. The correction involved the zero-point shift 
due to interstellar and atmospheric extinction and the slope 
due to the term associated to the second-order atmospheric 
extinction. At any given site with the same instrument and 
filter, this slope was assumed constant. Where relevant, we 
applied the correction. Remark that the effect is system- 
atically smaller in the F-band than in the B-band. The 
resulting differential magnitudes of every night were con- 
catenated per filter to produce the overall corrected B and 
V data sets. These data are presented in Tables |3] and SI 
The full tables are available in electronic form only. Table [T] 
also lists the numbers of observations collected at each ob- 
servatory. In the subsequent analysis, more than 7900 and 
2100 data points, respectively in the filters B and V , were 
treated. 

4. Simultaneous modelling of the light curves 

The B- and T^-data sets were used to plot the respective 
light curves phased against the orbital period of 1.7863748 
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Table 1. Campaigns and instrumentation. Total number of useful observations and hours. 



Year 


Site 


Country 


Tel. 


Period (& Obs.) 


Nights 


No. 


Hours 


Filter 


CODE 


2004 


Skinakas 


Crete 


1.3-m 


Jun-July (AS) 


4 


573 


17 


B 


SKBl 


2004 


Sobaeksan 


Korea 


0.61-m 


March (SLK) 


3 


379 


10 


B 


KRBl 


2005 


Skinakas 


Crete 


1.3-m 


May- Aug (AS) 


6 


1091 


29 


B 


SKB2 


2005 


Mt. Lemmon 


Arizona 


1.0-m 


Jun-July (SLK) 


5 


643 


21 


B 


KRB2 


2005 


Monegrillo 


Spain 


0.4-m 


Jun-July (JV) 


15 


1208 


67 


V 


SPVl 


2005 


Beersel 


Belgium 


0.4-m 


Jun-July (PVC) 


5 


430 


6 


B 


BHOB 


2005 


Beersel 


Belgium 


0.4-m 


Jun-July (PVC) 


2 


517 


6 


V 


BHOV 


2006 


Skinakas 


Crete 


1.3-m 


Jun-July (AS) 


9 


1427 


35 


B 


SKB3 


2006 


Sierra Nevada 


Spain 


1.5-m 


Aug (ER+MLG) 


9 


835 


18 


B 


SPBl 


2007 


Sierra Nevada 


Spain 


1.5-m 


Mar-May(ER-KMLG) 


8 


2158 


18 


B 


SPB2 


2007 


Beersel 


Belgium 


0.4-m 


March (PVC) 


1 


83 


3 


V 


BHOV 


2008 


NAO Rozhen 


Bulgaria 


2.0-m 


June- July (AS) 


2 


358 


12 


B 


ROZB 


2008 


Monegrillo 


Spain 


0.4-m 


June- July (JV) 


3 


372 


10 


V 


SPV2 


2004-08 


All 










2180 


86 


V 


ALL 


2004-08 


All 










7960 


166 


B 


ALL 



Table 2. Coordinates and standard system calibrated magnitudes with their errors for the target and comparison stars 



Ident 


CSC 


RA 


DEC 


B 


V 


R 


B-V 


V -R 


CT Her 


GSC01509-1142 


16 20 26.57 


+18 27 16.9 


11.33±0.03 


11.12±0.03 


11.00±0.03 


0.21±0.04 


0.12±0.04 


CI" 


GSC01509-1140 


16 20 33.89 


-hl8 27 18.5 


12.06±0.03 


10.80±0.03 


10.11±0.03 


1.26±0.04 


0.69±0.04 


C2''' 


GSC01509-1130 


16 20 18.61 


+18 27 32.6 


13.41±0.03 


12.40±0.03 


11.88±0.03 


1.01±0.04 


0.52±0.04 


C3 


GSC01509-0901 


16 20 43.00 


-hl8 30 53.2 


13.18±0.03 


12.67±0.03 


12.38±0.03 


0.51±0.04 


0.29±0.04 



" CI was also used by Kim et al. (|2004ap 



^ C2' is not the same star as C2 (Kim et al. I2004a|l 



Table 3. CT Her B light curve (first three lines, table 
available in electronic form) 



HJD 


B 


Code 


Year 


2453954.47391 


11.460156 


SPBl 


2006 


2453954.47529 


11.462120 


SPBl 


2006 


2453954.47609 


11.460165 


SPBl 


2006 



Table 4. CT Her V light curve (first three lines, table 
available in electronic form) 



HJD 


V 


Code 


Year 


2453557.38195 


11.953100 


SPVl 


2005 


2453557.38406 


11.938400 


SPVl 


2005 


2453557.38617 


11.921600 


SPVl 


2005 




0.4 0.6 
Orbital phase 



Fig. 1. Observed and theoretical light curve in the filter V . 
Observations are shown as circles. The solid line represents 
the model. 



days (jSamus et al. 2009)) . Using this ephemeris, we rede- 
termined the (single) epoch of primary eclipse observed in 
the filter V and obtained the improved value of HJDq — 
2442522.93270 ± 0.00006 days. The phased light curve in 
the filter V is almost complete (Fig. 1). However, we lack 
most of the primary minimum in the filter B (Fig. 2). Our 
current objective is to compute an adequate model from 
both curves in order to remove the geometric and photo- 
metric effects caused by the orbital motion and to extract 
the smallest possible residuals for a subsequent frequency- 
analysis. Concise reports of previous analyses (Lampens et 
al. I2008al I2008b|) did not include all of the above mentioned 
data sets. In particular, the full coverage of the primary 
eclipse was not (yet) achieved in the filter V. The data col- 
lected in 2008-2009 were thus mostly aimed at filling the 
gaps in orbital coverage. 



We performed a simultaneous modelling of the B and 

V light curves of CT Her using the light-curve fitting pro- 
gramme Phoebe, version 0.31a, in Mode 5, which corre- 
sponds to a semi-detached binary configuration in which 
the secondary component fills its limiting lobe, as required 
for the oEA stars. Phoebe (|Prsa &: Zwitter 2005bp is a 
package which enables to compute models of eclipsing bina- 
ries based on observed photometric and radial velocity data. 
It relies on the 2003-version of the widely used Wilson- 
Devinney code ( [Wilson &: Devinney 1971[ Wilson 119791 
I1990|) . First, we employed the almost complete and homo- 
geneous y-light curve obtained by one of us (JVS) during 
the years 2005 and 2008 to compute an initial model; next, 
using this initial model, we corrected the individual B and 

V data sets (where needed) for the differential effects of at- 
mospheric extinction as explained in Sect|21 Only then did 
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g 12.00 



12.50 



0.4 0.6 
Orbital phase 



Fig. 2. Observed and theoretical light curve in the filter B. 
Same legend as in Fig. 1. 
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Fig. 3. as a function of Ti and g 



we compute a best-fitting solution for the combined data 
series. 

Since CT Her is classified as an A3V star 
(jKuznetsova fc Svechnikov 1990p . we initially adopted 
8700 K as the surface temperature of the primary 
component. Except for the orbital ephemeris, all other 
parameters were set as adjustable parameters during the 
minimization procedure: this concerns the secondary star's 
surface temperature, T2, the inclination, i, the mass ratio, 
g, the dimensionless potential, fii, and the fractional 
luminosities of the primary component, Li, in B- and 
y- light (6 free parameters). We adopted the roughly 
determined absolute photometric and geometric elements 
(jSvechnikov fc Kuzne tsova 1990^ as starting values (e.g. 
i = 82° for the inclination and q = 0.27 for the mass ratio), 
together with an estimate of the secondary component's 
surface temperature, T2 = 5800 K, corresponding to a 
spectral type of G3IV. Other parameters such as the 
gravity darkening coefficients gi and 52 and the albedos 
Al and A2 were set to the theoretical values corresponding 
to a radiative atmosphere {gi — 1.0 and Ai — 1.0) in 
the case of the primary component and to a convective 
atmosphere (32 = 0.32 and A2 — 0.5, IRucmski 1969 ) in 
the case of the secondary component. The limb darkening 
coefficients in the B- and F-bands were taken from the 



tables by Van Hamme (1993) In all the runs, convergence 
was achieved after only two or three iterations using the 
differential correction method. 

By scanning the distribution of the function to be min- 
imized, x^, as a function of both Ti and the mass ratio. 




Orbital phase 



Fig. 4. Residual light curve phased against Porb (filter V). 



q, we found two regions of possible solutions: one region 
near Ti = 8700 K and another one near Ti = 8200 K. 
Fig. 3 illustrates the evolution of as a function of the 
parameters Ti and q. Thus, we considered two equivalent 
solutions in the subsequent analysis: one solution associ- 
ated with Ti = 8700 K (model 71) and one associated with 
Ti = 8200 K (model 72), with their corresponding values of 
T2 and q. The latter value of Ti corresponds to a spectral 
type of A5V and to a colour index {B — V)— 0.15 assuming 
zero reddening (cf. also the observed colour index {B — V) 
in Table[2])-The resulting parameter values and their formal 
uncertainties are listed in Table [5] Both models fit the ob- 
served B and V light curves very well, as evidenced by the 
small values of the x^ function and the 1-2%-level scatters 
of their residual data sets. The synthetic (and observed) 
light curves corresponding to the first model are illustrated 
by Figs. 1 and 2. The synthetic light curves corresponding 
to the alternative model are indistinguishable. 

The main difference between both models is the sur- 
face effective temperature of the secondary component, T2, 
which is shifted by about 200 K. The largest changes with 
respect to the initial values are found in the parameters 
T2 and q: both are smaller than their adopted first guesses 
(respectively 5800 K and 0.27). The mass ratio is signif- 
icantly smaller than previously assumed (i.e. q < 0.20). 
Nonetheless, even though the components' effective tem- 
peratures are different, both models are obviously similar 
since their derived physical parameters lie very close to 
one another (cf. the masses, gravities and radii of Table [TT|) . 

From here on, even though the values of the x^ function 
are slightly smaller in the former case (cf. Table [S]) , we will 
adopt the solution derived with Ti = 8200 K (model 72), 
the reason being that we obtained a better fit with this 
temperature to a few observed regions of the spectrum 
of CT Her (Sect. We will furthermore show that the 
frequency-analysis of the corresponding residuals is inde- 
pendent of the choice of either one of the two proposed 
models (Sect. [5]). 



5. Frequency analysis of the residual light curves 

The light curve model with Ti = 8200 K was subsequently 
subtracted from the original light curves of CT Her to 
search for short-period pulsations in the residual data. The 
phased residual light curves are plotted in Figs. 4 and 5, 
respectively in the filters V and B. The adjustment is 
excellent for the B-data set (without a primary minimum) 
and fair for the F-data set. Weights were associated based 
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Orbital phase 



Fig. 5. Residual light curve phased against Port (filter B). 



Table 5. Parameters of the simultaneous light curve solu- 
tions for CT Her, including their formal error. 



Parameter 


Model 71 
Filter B Filter V 


Model 72 
Filter B Filter V 




81.95 ± 0.01 


81.75 ± 0.01 


q 


0.1413 ± 0.0004 


0.1456 ± 0.0003 


Ti (K) 


8700" 


8200" 


T2 (K) 


4651 ± 7 


4489 ± 7 




4.406 ± 0.006 


4.426 ± 0.006 




2.079 


2.091 


[Li/{Li+L2)] 


0.969 0.932 


0.970 0.935 


[L2KU + L2) 


0.031 0.068 


0.030 0.065 


91 


1.0" 


1.0" 


92 


0.32" 


0.32" 




1.0" 


1.0" 


A2 


0.5" 


0.5" 




0.561'' 0.476'' 


0.585'' 0.505'' 




0.934'' 0.785'' 


0.962'' O.8II'' 


No. data 


7960 2180 
0.0160 0.0055 


7960 2180 
0.0163 0.0056 



" adopted value 

'' Van Hamme (1993) tables 

xi, X2 are the limb darkening coefficients. ^2 is the dimen- 
sionless potential of the secondary component. 



on the overall standard deviations of the (CT Her - CI) 
differential i?-magnitudes computed night by night: while 
1.0 was adopted in most cases, some sets with higher noise 
were allocated a relative weight of 0.2. This is the case for 
a few nights at the Beersel Hills Observatory (2005) and 
two nights at the Observatory of Sierra Nevada (2007). 
We restricted the residual data to the orbital phase bins 
between 0.05-0.95 in order to omit the phase of primary 
minimum (sparsely covered by our V-observations) . The 
fact that various outliers were found close to a primary 
minimum indicates that it is hard to model both effects 
simultaneously during this phase of rapid and steep 
light variation. On the other hand, the outliers are only 
associated to partially observed eclipses, whereas the fully 
observed primary eclipse of 2008 shows normal residuals. 
Such small effects might also have been introduced by the 
airmass-dependent corrections (since more prominent at 
the beginning/end of the night). After the removal of a few 
outliers larger than or equal to 0.02 mag, the remaining 
standard deviations are 6.4 mmag in the l^-band (with 



1958 residual data points) and 6.3 mmag in the B-band 
(with 7625 residual data points). 



5.1. Results from the B-data 

We performed Fourier analyses with Period04 
( |Lenz fc Breger 2005[ ). Fig. 6 illustrates the frequency 
search in the range 0-80 for the larger residual data 
set (_B-band): the spectral window is shown in the top 
panel while the following panels show the initial peri- 
odogram and the periodograms successively prewhitened 
of the strongest signal from each previous run. The last 
panel shows the residual periodogram with the mean noise 
level computed for frequency bins of width 5 d~^. We 
identified fifteen frequencies with an amplitude-to-noise 
ratio larger than or equal to 4.0 (the empirical criterion 
introduced by Breger et al. [Breger et al. 1993| . The 
frequencies, amplitudes, residual standard deviations, 
signal-to-noise ratios and the removed fraction of the 
initial variance, 1 — i? = {cfnit ~ '^res) / '^fnin obtained from 
a multi-parameter fit of the residuals to a solution with 
fifteen frequencies, are listed in Table [HI The errors on the 
frequencies and the amplitudes were computed using 120 
Monte Carlo simulations of Period04. The signal-to-noise 
ratios are based on the mean noise levels of the peri- 
odograms of the residuals. Remark that eight frequencies 
are concentrated in the range 43.5 - 53.5 d~^. The main 
frequency of 52.93664 d~^ has a semi-amplitude of 3.3 
mmag in B. This frequency corresponds to a periodicity 
of 27.2 min, and is the one found by Kim et al. (|2004a|l . A 
few oEA stars have even shorter main pulsation periods: 
e.g. RZ Cas (/i = 64.1935 d^^, Rodriguez et al.[2niH|) and 
AS Eri (/i = 59.03116 d"^ , Mkrtichian et al. l2004|) . 

All seven frequencies in the range / < 5 (i~^are caused 
by small imperfections linked to the reduction, extraction 
of the orbital variations and statistical fiuctuations in 
the data sets. Since they are irrelevant for our study, 
we kept them fixed during the Monte Carlo simulations. 
Some of the other eight frequencies found are affected 
by the 1 d~^ aliasing phenomenon: they are flagged in 
Table m At this point, we conclude that eight (pulsation) 
frequencies are found to be significant. After prewhitening 
for the multi-frequency solution, the remaining standard 
deviation in the _B-band equals 4.3 mmag, removing 
53% of the initial variance. High-quality residual i?-band 
light curves are presented in Fig. 7, demonstrating a good 
agreement between the residuals and the proposed solution. 



5.2. Results from the V-data 

Likewise, the complete set of the T^-band residual data was 
analysed. Since these data are less numerous and more 
noisy, less significant results were found. We identified ten 
frequencies with an amplitude-to-noise ratio larger than 
4.0, six of which are located in the range 43.5 - 53.5 d~^. 
The frequencies, amplitudes, residual standard deviations, 
signal-to-noise ratios and the removed fraction of the initial 
variance, 1 - R, resulting from a multi-parameter fit of the 
V^-band residual data set to a solution with ten frequencies, 
are listed in Table [71 All errors were computed using 120 
Monte Carlo simulations of Period04. Apart from the fre- 
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Fig. 6. 



20 40 60 

Frequency (cycles/day) 

Periodograms of CT Her computed during successive frequency analyses (filter B) 



quencies satisfying / < 5 d^^, the frequencies Gl = Fl, G3 
= F4, G6 = F8 and G9 — F14 were recovered, thus con- 
firming the presence of at least four of the eight previously 
detected frequencies. The dominant frequency has a semi- 
amplitude of 2.2 mmag and is located at 52.93762 d^^. It 
corresponds to the (1050 day)~^ alias frequency of Fl due 
to a 3-yr gap in the time series. With respect to the B- 
band analysis, two unidentified frequencies, G4 (at 49.497 
d~^) and GIO (at 48.068 d~^) were detected whereas the 
frequency F3 (at 47.600 d~^) was not confirmed. Whether 
these new frequencies might possibly be related to those al- 
ready found, remains unclear. Furthermore, the frequencies 
F6 and F15, which form an almost equally spaced triplet 
with G6 (= F8), were also not recovered. In Sect. 15. 3[ we 
will discuss these (apparently) partially incoherent results. 
After prewhitening for the solution with ten frequencies, 
the remaining standard deviation equals 5.1 mmag in the 
V-band. 



5.3. Variability of the pulsation amplitudes 

The full analyses in the B and y-bands were repeated af- 
ter subtraction of the alternative model using Ti = 8700 K 
(model 71). In the range above 24 d^^, the frequencies listed 
in Tables ini and [7] were confirmed, though with slight mod- 
ifications in the order of their appearance. This shows that 
the results of the Fourier analyses are perfectly independent 
of the adopted choice for the binary model ( owing to a) 
the fact that the binary models are almost identical and b) 
the different time scales involved). 

It would also seem that the results of the frequency 
analyses in both filters tend to a common solution, were it 
not for the presence of F3 (Table [H) and G4 (Table [H). In 
addition, the order in which the frequencies were detected 
in the respective data sets is not identical. Therefore, we 



Table 6. 

(model 72) 



_B-band frequency-analysis of the residuals 



ID 


Frequency 


Amplit. 


^ res 


S/N 1 - R 




(± error) 


(± 0.1) 










mmag 


mmag 




Fl 


52.93664 (± 0.8) 


3.3 


5.9 


12.7 


F2 


0.97682 (fixed) 


2.8 


5.6 


11.4 


F3 


47.59996 (± 2) 


1.8 


5.4 


7.2 


F4 


49.20822 (± 2) 


1.8 


5.3 


6.6 


F5 


2.37738 (fixed) 


2.0 


5.1 


8.1 


F6 


45.69130 (± 3)* 


1.3 


5.0 


6.2 


F7 


1.00032 (fixed) 


2.1 


4.9 


8.4 


F8 


45.44028 (± 2)* 


1.4 


4.8 


6.8 


F9 


3.44103 (fixed) 


1.1 


4.7 


4.3 


FIO 


53.23747 (± 2) 


1.2 


4.6 


4.8 


FU 


4.08614 (fixed) 


1.1 


4.5 


4.4 


F12 


0.61211 (fixed) 


1.2 


4.5 


4.7 


F13 


4.16902 (fixed) 


1.1 


4.4 


4.3 


F14 


43.76799 (± 3) 


1.1 


4.4 


5.4 


F15 


45.56395 (± 4)* 


1.0 


4.3 


4.6 0.53 



possibly affected by the 1 d ^aliasing 



investigated whether patterns could be found in the pul- 
sation amplitudes assuming that the frequency content is 
stable. All eight frequencies from Table [5] located in the 
range 43.5 - 53.5 d~^ were adopted to represent the fre- 
quency content of CT Her. Preference was given to the re- 
sults from the i3-band analysis because the aliasing is much 
stronger in the F-band data (particularly the 1 d~^ and the 
0.0018 d~^ aliasing effects). Next, we computed the multi- 
parameter solutions with the amplitudes and the phases as 
free parameters on a yearly basis. 
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-757.94 -757.90 -757.86 -757.83 -757.78 
Corrected JD 




-296.04 -296.00 -295.96 -395.92 -295.8 
Corrected J D 




-295.08 -295.04 -295,00 -294.96 -294.92 
Corrected JD 




47.50 47.54 47.58 47.62 47.66 47.70 47.74 47.78 
Corrected JD 





108.58 108.62 108.66 108.70 10S.74 

Corrected JD 



329.80 329.84 329.88 329.92 329.96 330.00 
Corrected JD 




330.80 330. 



330,88 330.92 330.96 331.00 
Corrected JD 




331,88 331,92 331.96 332,00 
Corrected JD 



Fig. 7. Sample figures of the B-hand residual data set and corresponding multi-frequency model. 



Table 7. V-band frequency- analysis of the residuals 
(model 72) 



Table 8, Semi- amplitudes of the pulsation frequenties {B- 
band) 



ID 


Frequency 
(± error) 


Amplit, 
(± 0,2) 




S/N 1 - R 






mmag 


mmag 




Gl 


52,93762 (± 5)*^ 


2,2 


6,4 


7,4 


G2 


0,56716 (fixed) 


2,2"= 


6,1 


6,4 


G3 


49,20636 (± 6)* 


1,8 


5,9 


6,8 


G4 


49,49661 (± 7) 


1,6 


5,6 


6,2 


G5 


3,68294 (fixed) 
45,43726 (± 479)* 


1,9 


5,5 


5,6 


G6 


1,1 


5,4 


5,3 


G7 


0,38478 (fixed) 


2,3 


5,3 


6,8 


G8 


2,08342 (fixed) 


1,6 


5,2 


4,7 


G9 


43,76795 (± 22) 
48,06813 (± 125)* 


1,1 


5,3 


5,4 


GIO 


1,1 


5,2 


4,4 0,36 



(possibly) affected by the 1 d ^aliasing 

: is the 0,00095 d"^ahas frequency of 5,93664 d"^(Fl, Tab, [6]) 
the 0,0018 d"^alias frequency of 49,20822 d"^(F4, Tab, [5]) 



ID 


Amp, 


Amp, 


Amp, 


Amp, 




2004 


2005 


2006 


2007-2008 




(± 0,3) 


(± 0,2) 


(± 0,1) 


(± 0.1) 




mmag 


mmag 


mmag 


mmag 


No, data 


951 


1704 


2262 


2511 


AT 


97 


61 


72 


466 


Fl 


3,0 


3,2 


3,4 


3,2 


F3 


2,2 (± 0,4) 


1,1 


2,2"= 


1,6'= 


F4 


2,4 


2,1 


1,3 


2,0^= 


F6 


2,2 


1,7 


0,7" 


1,7 


F8 


0,9 


0,6 


1,4"= 


2,2 


FIO 


1,8 


1,2 


1,2 


1,8 


F14 


0,2 


1,0 


1,5 


0,8 


F15 


1,3 (± 0,5) 


0,5 (± 0,5) 


1,2*= 


1,2= 



the error is ± 0,2 mmag 



the error is ± 0,3 mmag 



The results of the computations are shown in Table [8] 
and illustrated in Fig, 10, While the amplitude of the fre- 
quency Fl can be considered to be stable over the entire 
observing season, this is not true for all the frequencies. The 
frequencies F4 and F6 behave similarly (first decreasing in 
amplitude until 2006 and then increasing) , whereas the fre- 
quencies F3 and F8 show a maximum of their amplitude, 
respectively in 2006 and 2007, In contrast, the frequencies 
F14 and F6 show variability of their amplitude which oc- 



curs in anti-phase. This behaviour explains why the fre- 
quencies detected in Sect, I5,l land l5,2l are not identical: the 
detection of Gl = Fl, G3 = F4 and G4 = F8 may be un- 
derstood in the light of the 2005 (partially supplemented by 
the 2007-2008) i?-band semi-amplitudes. In conclusion, the 
amplitudes of these pulsation frequencies show evidence of 
variability on time scales of 1-2 years, perhaps even less, 
with the exception of the most dominant frequency which 
has a stable amplitude. 
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Fig. 8. Semi- amplitudes of the pulsation frequencies versus 
time (filter B) 

Table 9. Radial velocities of CT Her A (first three lines, 
table available in electronic form) 



HJD 


RV 


Error 


2454247.47535 


-44.21 


3.20 


2454247.49155 


-52.88 


2.21 


2454248.37028 


-1.60 


2.63 



6. Radial velocities 



High-resolution spectra of CT Her were gathered with the 
2-m telescope of the NAO, Rozhen, and the Coude spec- 
trograph equipped with an AT200 Photometries camera 
(with a SITe SI003AB chip and pixel size of 24 //m), and 
a spectral resolution of 0.19 A/pixel. Three regions were 
observed: H/3 and the regions around MgH (A = 448lA) 
and Fel (A = 5455 A). Forty-seven spectra were collected 
during 17 nights: 5 nights in May- July 2007, 5 nights in 
April- June 2008, as well as 7 nights in April- July 2009. 
Typical exposure times were 1200 and 1800 s. The aver- 
age signal-to-noise ratio S/N is w 30. The spectra were re- 
duced with standard IRAF procedures. The corresponding 
radial velocities were measured with the cross-correlation 
technique using synthe tic spectra obtained w ith the pro- 
gramme SPECTRUM dCray fc Corbally 2004D , and a grid 
of LTE-atmosphere models of solar-type chemical composi- 
tion (jCas telli & Kurucz 2003) (Table [5]) . Combined spectra 
having S/N > 100 were used to estimate the effective tem- 
perature of CT Her A: the best fit for two regions (H^, 
Mgll) was obtained with Tcff = 8250 K, log 5 = 3.7 and 
usini = 60 km/s (Fig.lH]). 

Eight spectra were also collected with the FOCES spec- 
trograph attached to the DSAZ 2. 2-m telescope of the Calar 
Alto Observatory (CLA) in 2009. The resolution is 40600. 
The spectral region covers 5000-6000 A. The average S/N 
is « 25. The radial velocities of both components were 
measured with the cross-correlation technique using syn- 
thetic spectra generated for temperatures of 8200 K with 
log g — 4.0 and 4500 K with log g — 3.5, and for 
usini = 50 km/s. In this case, the 2-D cross-correlation 
programme TODCOR was used (jMazeh fc Zucker 1994|) . All 
radial velocity measurements are presented in Table HI 
available in electronic form only. We first checked that the 
radial velocities obtained for component A were compati- 
ble with the NAO radial velocities. In a plot versus time 
(cf. top panel of Fig. fTU]) . we detected a long-term variation 
which had to be removed before fitting. 



^ 1.0 

1 0.8 
c 



Ugll 4481 Region 



S/N=146 



■ T^i = 8750 I I 



S/N=173 






_T^, = 8750\,X\ 
-T^i = 8250 






4550 4600 4800 

Wovelength [A] 



Fig. 9. Top panel: combined observational spectra and syn- 
thetic spectra showing model A (Toff = 8250 K), and 
model B (Teff = 8750 K). Other panels: difference spectra 
in the sense observed minus model A (middle) and observed 
minus model B (bottom). 



A Fourier analysis of the residuals obtained from the 
raw velocities performed with Period04, and a prelimi- 
nary fit with velocity semi-amplitude Ka — 23.4 km/s and 
systemic velocity 7 ~ —15.8 km/s indicated a periodic- 
ity of K. 125.3 ± 2.2 days with semi-amplitude Kc — 11.4 
km/s. The corresponding ephemeris can be described by: 
HJD^ax = 2454157.6 + 125.3 x E. Fig. [TO] shows the 
model (solid line in top panel) and the raw and corrected 
radial velocity curves (bottom panels). The standard devia- 
tions of the radial velocity residuals amount to 11.69 km/s 
and 7.63 km/s, respectively for the raw data and for the 
data after subtraction of the model. 

The corrected radial velocity curves were next used 
in combination with the B and V light curves to en- 
able fitting of the following parameters with Phoebe: the 
semi-axis major expressed in solar radii. A, and the sys- 
temic velocity, 7. The new parameter values and their 
formal uncertainties are listed in Table [TOl This value 
of the semi-axis major, A, lies close to its approximate 
value (jSvechnikov &: Kuznetsova 1990p . The values of the 

function are slightly worse than before. This new solu- 
tion fits both light curves a little worse, particularly the 
primary eclipse in the V-filter, but accomodates the radial 
velocities reasonably: the means with their standard devi- 
ation of the residual radial velocities equals 1.7 ± 5.8 km/s 
and 1.4 ± 12.3 km/s, for component A and B respectively 
(Fig. [TT|1 . The main differences with the former solutions 
are the smaller value of the mass ratio (g = 0.127) and 
the higher value of the inclination [i — 82.8°). Table [TT] 
summarizes the physical properties of both components re- 
sulting from the combined data fitting. These properties 
appear to be compatible with the status of component A 
as a pulsator of type 5 Scuti. The computed masses are 
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Table 10. Parameters of the combined light and radial 
velocity solution for CT Her. 



Long-term RV curve of CT Her 



Parameter 



Filter B 



Filter V 



in 




82.75 ± 0.04 


q 




0.1267 ± 0.0006 


A 




8.48 ± 0.03 


7 {km/u 




-16.2 ± 0.2 


Ti (K) 




8200" 


T2 (K) 




4468 ± 7 


^1 




4.269 ± 0.008 


^2 " 




2.038 




+ L2)] 


0.969 0.932 


[L2/{Li 


+ L2)] 


0.031 0.068 


gi 




1.0" 






0.32" 






1.0" 


A2 




0.5" 


xi 




0.585* O.SOS** 


X2 " 




0.962'' O.8II'' 



No. data 
No. RV data 



7960 2180 

0.0162 0.0055 

55 (comp A) 7 (comp B) 

0.0038 0.0151 



" adopted value 

'' Van Hamme (1993) tables 

xi, X2 are the limb darkening coefficients. is the dimen- 
sionless potential of the secondary component. 



Table 11. Tentative absolute parameters of CT Her 



Parameter 


Comp A 


Comp B 


Mass (Mq) 


2.28 ±0.01 


0.29 


±0.04 


Radius (R©) 


2.06 ±0.06 


1.87 


±0.08 




8200 


4468 ± 80 


log g 


4.17 ± 0.02 


3.36 


± 0.02 


Luminosity (L©) 


17.4 ±2.4 


1.2 


±0.2 


Mboi (mag) 


1.65 ±0.15 


4.50 


±0.21 



furthermore consistent with those preliminarily derived by 
Hoffman fc Harrison (2009) from spectroscopy only (they 
obtained Mi 



2.31 ±0.11 Mq and M2 = 0.31 ± 0.09 M©, 
with g = 0.13 (in agreement with the q- values of Tables [5] 
and [TU|) ) . We propose that this new solution illustrates the 
uncertainty still remaining in the characterization of the 
physical parameters of CT Her: more realistic errors would 
thus correspond to the formal errors of Table [S] (Table [TU)) 
multiplied with a factor of 25-50 while an uncertainty of 
about 80 K would be closer to the truth in the case of 
T2. The data collected at the phase of primary minimum 
are probably affected by the pulsations (this is indeed the 
case for other oEA stars), which may introduce some in- 
determination in the light curve modelling. However, our 
radial velocity data, in particular of component B, are too 
scarce to enable an accurate determination of the mass 
ratio. Accurate component radial velocities obtained from 
high-resolution echelle spectra would be necessary to obtain 
a more consistent determination of the absolute parameters 
of CT Her. Meanwhile, the solution(s) derived in Sect.|3]are 
currently the most adequate one(s). 
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Fig. 10. Radial velocities of CT Her A obtained from NAO 
spectra (plusses) and CLA spectra (diamonds): a 125.3 days 
periodicity was removed before fitting in combination with 
the B- and V^-light curves. 



O Comp A 
• Comp B 



Orbital phase 



Fig. 11. Residual component radial velocities of CT Her. 



7. Possible effect on the pulsation analysis 

The presence of a third body can cause a cyclic variation of 
the systemic velocity. In that case, the resulting light travel 
time (LTT) effect will introduce periodical time delays in 
the pulsational analysis. We examined if any phase shifts 
could be detected by subdividing the S-filter data into 10 
subsets arranged according to the orbital phase bin asso- 
ciated with the long-term ephemeris discussed in Sect [S) 
We used the multivariate analysis method implemented in 
PERIOD04 dBreger 20051 ) to recompute the best-fitting 
phases for each subset with respect to the proposed multi- 
frequency solution. The outcome of these computations is 
that we find no obvious shift with respect to the initially 
adopted phases. Fig. [H] represents the computed phase 
shifts, and illustrates that there is no detection of a sys- 
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while the squares show the 0-C values using the ephemeris 
corresponding to the best-fit parabola: 



> Phase shift Fl 
. Phase shift F3 

> Phase shift F4 



LTT phase bin (orb. period ~ 125 days) 

Fig. 12. Phase shifts computed for the pulsation frequen- 
cies Fl, F3 and F4 as a function of the orbital phase bin 
(with respect to the long periodicity). 

Table 12. 0-C values for two models representing the 
evolution of Porh (first three used lines, table available in 
electronic form) 



HJD 


E 


O-Ci 


O-C2 


p/s 


Tech. 


2438894.7540 


-2031 


-0.0508 


-0.0337 


P 


vis 


2438953.7280 


-1998 


-0.0271 


-0.0103 


P 


vis 


2439248.5090 


-1833 


0.0020 


0.0175 


P 


vis 



tematical periodical pattern in the phases of the frequencies 
Fl, F3 and F4. It is clear that a lot of data in each subset 
are necessary to be able to reliably use the result of the 
multi-frequency fit. That is why we considered 6 out of 10 
subsets only. For these frequencies (Fl, F3 and F4) and 
for these subsets, the errors on the phases are sufhciently 
precise. 

Alternatively, one can derive a model making allowance 
for a periodic variation of the phase shifts of the multi- 
frequency solution. This approach has been applied as well 
and the results confirm that no significant differences on 
the (pulsation) frequencies are detected even though a 
slightly better overall fit was obtained, corresponding to 
a 129.6 days LTT (solution with Cres — 4.2 mmag). 



8. The 0-C diagram 

The same LTT effect also implies a cyclic variation of 
the orbital period. Therefore, we investigated the 0-C di- 
agram of CT Her. We used the list of observed times 
of light minima from the database of the Variable Star 
and Exoplanet Section of the Czech Astronomical Society 
(|2009p , from which we selected 146 observed times collected 
between 1965 and 2009 {JD > 2438761), after elimina- 
tion of the 26 first epochs. Then, we added two of our 
own observed minima (one minimum was not useful for 
the modelling of the B-band light curve because of an un- 
explained sudden drop during the eclipse and another iso- 
lated minimum was obtained by D. Litvinenko in an un- 
reported filter): 2453189.3784 ± 0.0002 {E = 5971) and 
24535 14.4961 ±0.0001 {E = 6153). Using Porb = 1.7863748 
days (jSamus et al. 2005)) . we recomputed the 0-C values 
and searched for the best representation using either a lin- 
ear or a parabolic model. Fig. [T51 illustrates both cases: the 
circles show the 0-C values with the known orbital period 



HJD„ 



-- 2442522.92914+ x 1^7863799 
^E^ X 7?394 X 10"^°. 



Such a model implies a larger orbital period with a de- 
crease dP/dt of 8.28 X 10^-^° days d~^{i.e. a secular pe- 
riod decrease dP/(Pdt) of 169 x lO^^yr"^). The rms is 
0.01049 days (906 s). To check its significance with respect 
to the linear model, we used an upper one-tailed F-test to 
verify whether the null hypothesis that the sums of squared 
residuals of both models are equal holds. The test statistic 
equals 



iix. 



2 

lin 



Xlar) X (n - 3))/{xlar X (3 - 2)) = 23.44. 



Since this value is much larger than the tabulated value 
Fo.oi(l, 145) = 6.81 fPezz uUo 20"TT]) . we reject the null hy- 
pothesis at the 1% significance level (type I error), mean- 
ing that the parabolic model is stati stically significant. We 
remark that [Soydugan et al. (2008) used all previous 0-C 
data (including the 26 first measurements) to derive an 
overall increase with a cyclic variation for the orbital period 
of CT Her. For evident reasons, we prefer to discuss only 
the well-observed part of this diagram. 

Table [T^ lists the cycle numbers and the associated 
residual values for both models. However, the first reported 
times do not fit these models at all as is also obvious from 
the 0-C diagram in the Atlas of 0-C diagrams of Eclipsing 
Binary Stars (|Kreiner et al. 20M]) . Moreover, the scatter 
remains large in both cases. Interestingly, if we restrict 
the observations only to the 24 most recent times observed 
using the CCD technique, we obtain a tight-fitting linear 
model with a high correlation coefficient (0.9) indicating an 
increased orbital period of 1.7863789 days with respect to 
the known period (cf. the speckled circles in Fig. [T3|) . The 
rms in this case is only 0.00135 days (117 s). 

Using the information about a probable variation in 
systemic velocity with a period of 125.3 days and a semi- 
amplitude of 11.4 km/s (Sect.H]), we estimate a mass func- 
tion equal to 



f{M) = 1.0385 X 10"^ X (1 - e^f/^ x K^P 



0.019 Mp 







and a corresponding 0-C amplitude of 131 s. Because this is 
of the same order as the computed rms (117 s), we conclude 
that the LTT effect is barely observable from the currently 
available data and that the 0-C diagram of CT Her does 
not contradict the proposed model of a variable systemic 
velocity due to the presence of a third body. Remark that 
the first two epochs of light minimum observed with a CCD 
show a fiat trend, in full agreement with the period men- 
tioned in the literature. 



9. Summary and conclusions 

CT Her is yet another oEA star, after Y Cam, in which a 
large set of pulsation frequencies has been identified thanks 
to a long-term photometric monitoring performed in two 
filters We confirm the presence of rapid pulsations with 
a dominant frequency of 52.93664 d~^. As many as eight 
significant frequencies, all of which are located in the range 
43.5-53.5 d~^, appear to reproduce the observed, complex 
behaviour reasonably well throughout the years 2004-2008. 
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o P= 1.7863748 d 

— larger P incl dP/dt 

^ P^ 1.7863800 d + dP/di 
linear fit after siibtracdon 

• CCD only 
larger P (CCD only) 



3500 5000 
Number of cycles elapsed 



tween pulsation and mass transfer in close binary systems, 
including the effects of tidal interaction. 
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Fig. 13. Evolution of the orbital period: 0-C diagrams of 
CT Her using a constant and a variable model. 



Furthermore, except for the most dominant frequency, 
we showed that the amplitudes of these pulsation frequen- 
cies are variable on time scales of 1-2 years, perhaps even 
less than one year, assuming that the frequency content 
is stable (the fact that several frequencies are commonly 
found in different data sets suggests this). Considering the 
absolute parameters of the (actual) primary component of 
CT Her (Table [TT|) . we conclude that the dominant fre- 
quency cannot correspond to the fundamental radial mode: 
rather, such a short period of pulsation for a 2.5 M© star is 
typical of a high overtone radial mode or a non-radial pul- 
sation mode of high order "p" (similar to RZ Cas, AS Eri 
or TW Dra). Using the improved ephemeris, we derived an 
accurate value of 94.56 for the ratio Porb/^puis- 

The analysis of recent spectra providing complementary 
radial velocities suggests that the light and radial velocity 
curves can be modelled simultaneously, provided that the 
radial velocity measurements are corrected for a long pe- 
riodicity of w 125 days with a semi-amplitude Kq = 11.4 
km/s. The presence of a third body in the system could ex- 
plain such a light travel time effect. However, we were not 
able to confirm this suggestion from an investigation of the 
0-C diagram of CT Her as the expected 0-C amplitude 
is of the same order as that of the noise. Neither did we 
detect an obvious orbital modulation in the phases of the 
multi-frequency solution. 

Although the match between the model and the B- 
and y-light curves is excellent, we cannot yet claim to 
know the system parameters of CT Her uniquely, especially 
because accurate component radial velocities are lacking. 
Such component radial velocities are much needed in or- 
der to obtain a consistent determination of the absolute 
parameters of this Algol-type binary. Meanwhile, we ac- 
quired a new series of high-resolution spectra of CT Her 
with the Hermes spectrograph attached to the Mercator 
telescope (jRaskin et al. 20lT|l . High signal-to-noise spectra 
should also be very useful for application of the technique o f 
spectra disentangling (e.g. TW Dra, Lehm ann et al. 2009^ . 
allowing to recover the component spectra from the com- 
posite ones and to study the residual line profile variations 
for a spectroscopic identification of the excited modes. The 
results obtained from this study show that CT Her is a 
particularly interesting binary system. oEA stars indeed 
provide extremely useful asteroseismic targets since they 
offer the opportunity to study in detail the connection be- 
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